Electrical conduction processes in ZnO in a wide temperature range 20—500 K 
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We have investigated the electrical conduction processes in as-grown and thermally cycled ZnO 
single crystal as well as as-grown ZnO polycrystalline films over the wide temperature range 20- 
500 K. In the case of ZnO single crystal between 110 and 500 K, two types of thermal activation 
conduction processes are observed. This is explained in terms of the existence of both shallow 
donors and intermediately deep donors which are consecutively excited to the conduction band 
as the temperature increases. By measuring the resistivity p(T) of a given single crystal after 
repeated thermal cycling in vacuum, we demonstrate that oxygen vacancies play an important role 
in governing the shallow donor concentrations but leave the activation energy (~27±2 meV) largely 
intact. In the case of polycrystalline films, two types of thermal activation conduction processes 
are also observed between ~150 and 500 K. Below ~150 K, we found an additional conduction 
process due to the nearest-neighbor-hopping conduction mechanism which takes place in the shallow 
impurity band. As the temperatures further decreases below ~80 K, a crossover to the Mott variable- 
range-hopping conduction process is observed. Taken together with our previous measurements on 
p(T) of ZnO polycrystalline films in the temperature range 2-100 K [Y. L. Huang et al, J. Appl. 
Phys. 107, 063715 (2010)], this work establishes a quite complete picture of the overall electrical 
conduction mechanisms in the ZnO material from liquid-helium temperatures up to 500 K. 

PACS numbers: 72. 20. My, 71.30.Tm, 73.50.Jt 



I. INTRODUCTION 

Zinc oxide (ZnO) is a wide band gap semiconductor 
with a direct energy gap of ss3.4 eV at room tempera- 
ture. This class of material is usually unintentionally (or, 
natively) doped with n-typc impurities, such as oxygen 
vacancies, Zn interstitials, and impurity hydrogen atoms. 
However, the exact doping behavior of different types of 
unintentional and intentional (e.g., group I or group III) 
dopants is still under much theoretical and experimental 
debate." 

ZnO has recently attracted intense attention owing to 
its potential prospects in optoelectronic applications^— 
Compared with the widely studied optical properties, the 
electrical conduction processes in ZnO materials (single 
crystals, polycrystalline films, etc.) have not much been 
investigated over a wide range of temperature^— Phys- 
ically, the measurements on the temperature dependence 
of resistivity, p(T), can provide very useful information 
on the carrier transport mechanisms as well as the asso- 
ciated impurity levels and their distributions in energy 
in a given semiconducting material. Recently, Huang 
et a/.— have investigated the electrical conduction pro- 
cesses in a series of oxygen deficient polycrystalline ZnO 
films (~1 pm thick) in the temperature range 2-100 
K. They observed the three-dimensional Mott variable- 
range-hopping (VRH) conduction process below ^100 K. 
As the temperature further decreased below ^25 K, a 
crossover to the Efros-Shklovskii VRH conduction pro- 
cess was found. These two types of VRH conduction 
processes in ZnO films have been successfully explained 



in a coherent manner. On the other hand, the charge 
transport properties above 100 K has not been explic- 
itly addressed in that work, due to the lacking of p(T) 
data above 300 K. In this work, we aim at extending 
the measurements on p(T) to temperatures up to 500 
K in order to further unravel the overall electrical con- 
duction mechanisms in ZnO, especially at T>100 K. We 
have measured an as-grown ZnO single crystal and re- 
measured the series of as-grown ZnO polycrystalline films 
which were first studied in Ref. [Til The as-grown single 
crystal has been repeatedly thermally cycled in vacuum 
and their p(T) has been measured during each thermal 
cycling process. The polycrystalline films had been de- 
posited in different oxygen atmospheres to contain differ- 
ent amounts of oxygen vacancies so that they possessed 
differing shallow donor concentrations when as prepared. 
Our new results are reported below. 

We would like to emphasize that all the p(T) mea- 
surements reported in this study have been carried out 
in a vacuum (<10 -4 torr) and in a "dark" environment 
without involving any photo-induced carriers. There- 
fore, complex thermal annealing effects which might re- 
sult from resistivity measurements performed in, e.g., 
an air (or an active gas) atmosphere at highly elevated 
temperatures,— can be safely avoided. 

This paper is organized as follows. In Sec. II, we briefly 
discuss our experimental method for sample fabrication 
and resistance measurements. Our results of the temper- 
ature dependence of resistance in the wide temperature 
interval 20-500 K and their physical interpretations are 
presented in Sec. III. Section IV contains our conclusion 
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as well as a summary of the overall electrical conduction 
mechanisms in the ZnO materials. 



II. EXPERIMENTAL METHOD 

Our ZnO single crystal (2.52x1.94x0.45 mm 3 ) were 
grown by the hydrothcrmal method and was obtained 
from a commercial supplier ^ Unfortunately, the detailed 
growth conditions were unavailable. However, the p{T) 
data do provide us meaningful information on the car- 
rier transport processes in the ZnO material (see below) . 
This as-grown crystal has a room temperature resistiv- 
ity of «5.45xl0 4 fi cm. Our ZnO polycrystalline films 
were fabricated by the standard rf sputtering deposition 
method on glass substrates. The starting ZnO target was 
of 99.99% purity, and the glass substrates were held at 
550° C during the deposition process. The films were de- 
posited in varying mixtures of argon and oxygen gases. 
The p(T) of these films had previously been measured 
and discussed by Huang et al^ for T<100 K. 

In this work, four-probe resistance measurements were 
carried out using a Keithley Model K-6430 as a current 
source and a high- impedance (Tf2) Keithley Model K-617 
as a voltmeter. The samples were mounted on the sam- 
ple holder which was situated inside a stainless vacuum 
can in a JANIS CCS350 closed-cycle refrigerator (10- 
500 K). The working pressure in the vacuum can was 
<10 -4 torr throughout all resistivity measurement pro- 
cesses, whether the vacuum can was connected to a diffu- 
sion pump (at T>250 K) or not (at T<250 K, see below 
for further discussion). Moreover, the stainless vacuum 
can contained no optical windows, and thus our measure- 
ments were performed in the dark, i.e., this experiment 
did not invoke any kind of photo-induced currents. The 
four-terminal current leads and voltage leads (pogo pins) 
were attached to the samples with the aid of silver paste. 
It should be noted that the resistances reported in this 
work were all measured by scanning the current-voltage 
(I- V) characteristics at various fixed temperatures be- 
tween 20 and 500 K. The resistance at a given temper- 
ature was determined from the regime around the zero 
bias voltage, where the I-V curve was linear. In fact, our 
I- V curves in every sample were linear over a wide range 
of bias voltage. The inset of Fig. [T] shows the linear I-V 
characteristics of the as-grown ZnO single crystal at four 
representative temperatures, as indicated in the caption 
to Fig. [TJ 



III. RESULTS AND DISCUSSION 

A. ZnO single crystal: thermal activation 
conduction over 110—500 K 

Figure [1] plots the measured resistivity as a function of 
reciprocal temperature for the as-grown single crystal be- 
tween 110 and 500 K. (The resistance below ~110 K be- 



came too large to be accurately measured.) Two distinct 
slopes are explicitly seen in the two different tempera- 
ture regimes of ~110~220 K (straight dotted line) and 
^290-500 K (straight dashed line), respectively. Such 
Arrhenius-type behavior immediately suggests that the 
responsible charge transport arises from the thermal ac- 
tivation conduction processes. Indeed, we found that our 
measured p(T) data can be quantitatively described by 
the equation 

p(T)- 1 = p- x l e- E ^' kBT + p^ e - E -' kBT , (1) 

where p\ and pi are temperature insensitive resistivity 
prefactors, and E\ and £2 are the relevant activation en- 
ergies associated with the two kinds of thermal activation 
conduction processes, is the Boltzmann constant. In 
Fig. 1, the solid curve is a least-square fit to Eq. (p}. It 
clearly manifests that Eq. (p} can well describe the ex- 
perimental result. Our fitted values of p\. E\, P2, and 
E2 are listed in Table I. This observation illustrates that 
there exist in the as-grown single crystal a group of shal- 
low donors with an activation energy of E2~29 meV and 
a group of intermediately deep donors with an activation 
energy of i?i«330 meV. (In this work, we shall maintain 
the convention E\ >E^ .) At the sufficiently high measure- 
ment temperatures of T>290 K, both the shallow and 
the intermediately deep donors are excited to the con- 
duction band and are responsible for the electrical trans- 
port. However, the exponential temperature dependence 
of p is mainly governed by the number of those interme- 
diately deep donors which are being excited to the con- 
duction band (the so-called ^-conduction channel). As 
T reduces below ^220 K, the intermediately deep donors 
become essentially intact while the shallow donors can 
still be readily excited to the conduction band. The T 
dependence of p is thus largely determined by the number 
of those shallow donors being excited to the conduction 
band (the so-called ^-conduction channel). In short, 
Eq. (fTJ) represents the "band conduction" processes. 

It should be noted that our extracted E\ and E2 values 
are in good accord with the corresponding activation en- 
ergies previously obtained by Look^. (-~340 meV), Wenck- 
stern et alr^ (~34-37 and -~300-370 meV), and Schifano 
et al*£- (~30, ^50, and ~290 meV) in various ZnO single 
crystals and in different temperature intervals. Recently, 
in a series of single-crystalline ZnO nanowires, Chiu et 
alJ£ (Tsai et alM-) have also reported an E2 activation 
energy of ~30-40 (~25) meV. In other words, we may 
conclude that different ZnO single crystals grown by dif- 
fering methods and under differing conditions possess es- 
sentially similar types of shallow and intermediately deep 
donor levels. (Obviously, the donor concentrations would 
vary from experiment to experiment.) 

In Fig. [5J we plot our measured resistivities as a func- 
tion of reciprocal temperature for the as-grown and then 
repeatedly thermally cycled single crystal. The closed 
squares (taken from Fig. [1]), open squares, closed cir- 
cles, and open circles represent the as-grown, first-cycled, 
second-cycled, and third-cycled samples, respectively. In 
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TABLE I. Relevant parameters of as-grown and thermally cycled ZnO single crystal, pi, Ei, p 2 , and E2 are defined in Eq. (1). 



single crystal 


p(300K) (fi cm) 


pi (fi cm) 


E ± (meV) 


p2 (SI cm) 


E 2 (meV) 


as-grown 


5.45xl0 4 


0.144 


334 


1.47x10 s 


29.2 


first-cycled 


4.56xl0 4 


0.154 


330 


2.50xl0 4 


28.1 


second-cycled 


2.70xl0 4 


0.145 


332 


1.72xl0 4 


24.7 


third-cycled 


2.01xl0 4 


0.125 


337 


1.02xl0 4 


26.4 




1/7X1/K) l/r(l/K) 



FIG. 1. (color online) Logarithm of resistivity as a function 
of reciprocal temperature for as-grown ZnO single crystal be- 
tween 110 and 500 K. The symbols are the experimental data 
and the solid curve is a least-squares fit to Eq. JTJ). The 
straight dashed and dotted lines are guides to the eye. Inset: 
Current-voltage curves at four temperatures of 490, 410, 370, 
and 130 K. Notice that the I-V characteristics are linear, i.e., 
ohmic. 



practice, our p(T) data for every sample were measured 
from 500 K down. The sample was placed in the dark 
vacuum can and first heated to 500 K. The T dependence 
of p was then measured with progressive lowing of tem- 
perature. During the data acquisition process from 500 
K down to 250 K, the dark vacuum can was connected to 
a diffusion pump for continuous pumping. The diffusion 
pump was disconnected when the temperature reached 
250 KJ£ Such a measurement process (which lasted for 
^270 min between 250 and 500 K) had effectively served 
as a "thermal annealing" (in vacuum) process. There- 
fore, the oxygen contain in the sample was reduced from 
that in the previous run. In other words, the amount 
of oxygen vacancies increased, leading to a decrease in 
the sample resistivity. We have repeated the p{T) mea- 
surements in the temperature interval 110-500 K for four 
times, and denoted the sample by as-grown, first-cycled, 
second-cycled, and third-cycled samples, as mentioned. 

Inspection of Fig. [2] indicates that, in the temperature 
interval ~330-500 K, the p(T) data collapse closely onto 



FIG. 2. (color online) Logarithm of resistivity as a function 
of reciprocal temperature for the as-grown (closed squares, 
taken from Fig. [lj, first-cycled (open squares), second-cycled 
(closed circles), and third-cycled (open circles) ZnO single 
crystal. The solid curves are least-squares fits to Eq. fTJ, see 
text. 



a straight line. On the other hand, below ~330 K, the 
measured p(T) data systematically decrease after each 
repeated thermal cycling. Quantitatively, the resistivity 
ratios of the as-grown sample to the third-cycled sample 
are ~1.0 at 460 K, ~2.7 at 300 K, and ~11 at 220 K. This 
variation in resistivity ratio suggests that the rich amount 
of the intermediately deep donors^ and their associated 
activation energies, i.e., the E\ values, are barely affected 
by our measuring and/or thermal cycling process. This 
is expected, since the activation energies associated with 
these donor levels are comparatively large while our ther- 
mal cycling temperature is relatively low. 

What is more interesting is the T behavior of p between 
110 and 220 K. While the p value at a given T decreases 
with repeated thermal cycling as mentioned, the temper- 
ature dependence of p maintains essentially unchanged. 
Numerically, our least-squares fitted slopes, i.e., the £2 
values, in this temperature interval vary only by ~10%, 
see Table I. This quantitative result is meaningful. It 
strongly implies that the amount of oxygen vacancies, 
namely, the shallow donor concentration, in the sample 
is increased after each repeated thermal cycling. Never- 
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theless, the thermal activation energy, namely, the shal- 
low donor level below the conduction band minimum, 
remains essentially unchanged. (See Ref. for a brief 
discussion on the thermal-cycling-induced weak disper- 
sion of the shallow donor energy levels.) This observation 
is strongly supportive of the important roles of oxygen 
vacancies as the shallow donors in ZnO. On the theoret- 
ical side, based on the formation energy consideration, 
it has recently been argued that oxygen vacancies would 
form negative-U centers, rather than shallow donors, in 
the ZnO material £lr— This theoretical prediction is not 
supported by the present result, Fig. [5J This puzzling 
issue requires further clarification. 



B. ZnO polycrystalline films 

1. Thermal activation conduction and 
nearest-neighbor-hopping conduction: 90-500 K 

Figure [3] shows the variation in the logarithm of re- 
sistivity with reciprocal temperature for four ZnO poly- 
crystalline films between 90 and 500 K, as indicated. In 
sharp contrast to the case of single crystal (Fig. [5]), Fig. [3] 
clearly reveals that p smoothly increases with decreas- 
ing temperature. That is, there does not exist a visible 
straight regime in any temperature interval in the logp- 
(1/T) plot from 500 down to 90 K. Quantitatively, we 
found that these results can not be described by Eq. (fTJ . 
Instead, they can be described by the following equation 

p(T)- 1 = p-l e -Ei/k B T + p -l e -E 2 /k B T + p -l e -E 3 /k B T ^ 

(2) 

where pi, E\, p2, and E 2 have the similar meaning as 
defined in Eq. (TTJ) . The third term on the right hand side 
represents a new, additional electrical conduction chan- 
nel which is characterized by a temperature insensitive 
resistivity prefactor p% and a thermal activation energy 
E3. (In the following discussion, we shall maintain the 
convention Ei>E 2 >E 3 .) 

Our fitted results with Eq. ^ are plotted as the solid 
curves in Fig. [3] and the values of the adjustable pa- 
rameters are listed in Table II. We obtain the thermal 
activation energies -Ei«113-135 meV, i?2~26-39 meV, 
and i?3~2.3-5.7 meV. The three types of conduction pro- 
cesses dominate the T dependence of p in the tempera- 
ture intervals of ~400-500 K (the E\ -conduction chan- 
nel), ^150-400 K (the ^-conduction channel), and ~90- 
150 K (the ^-conduction channel). Physically, the inl- 
and the ^-conduction processes are just those described 
in the above subsection. However, the E\ activation en- 
ergy is markedly reduced from «330 meV for single crys- 
tal to «125±10 mev for polycrystalline films (see further 
discussion below) . The additional ^-conduction channel 
can be ascribed to the nearest-neighbor-hopping (NNH) 
conduction mechanism which takes place in the shallow 
"impurity band." That is, as the temperature decreases 
below ~150 K, even the shallow donors (not mentioning 
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FIG. 3. Logarithm of resistivity as a function of reciprocal 
temperature for four ZnO polycrystalline films between 90 
and 500 K, as indicated. The symbols are the experimental 
data and the solid curves are least-squares fits to Eq. 



the intermediately deep donors) can no longer be effec- 
tively excited to the conduction band. The charge trans- 
port is then largely governed by those electrons which 
hop from an occupied state (which lies below the Fermi 
energy) to an unoccupied nearest-neighbor state (which 
lies above the Fermi energy) in the impurity band. Such 
a hopping process is phonon assisted and requires a small 
amount of energy, i.e., the E3 thermal activation energy. 
It should be noted that a shallow impurity band can form 
in the oxygen deficient ZnO polycrystalline films but not 
in the ZnO single crystal, because the former samples 
contain far higher amounts of oxygen vacancies. A large 
amount of oxygen vacancies causes a notable dispersion of 
the donor levels or ionization energies, and thus forming 
an impurity band.— Due to the slight self-compensation 
characteristic of unintentionally doped ZnO; 16 ' 26 the im- 
purity band is partially filled, making the NNH conduc- 
tion process possible. Based on this picture, the £3 value 
should decrease with increasing oxygen vacancies, be- 
cause there would be more randomly distributed donor 
levels per unit energy per unit volume. Indeed, inspec- 
tion of Table II reveals that our extracted E 3 value (on 
average) decreases with the decreasing O2 flux applied 
during the fabrication process. 

In contrast, in the case of single crystal, the shallow 
donor concentration is sufficiently low and thus, on aver- 
age, the donor impurities should lie sufficiently far apart 
in space and in energy. Therefore, the NNH conduction 
process can only play a negligible role, as compared to 
the £"2 -conduction process, even at temperatures down 
to 110 K. Inspection of Tables I and II indicates that the 
p(300 K) values of our polycrystalline films are more than 
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TABLE II. Values of relevant parameters for six oxygen deficient ZnO polycrystalline films, pi (i=l, 2, 3) and Ei are defined 
in Eq. (2). Notice that these samples are taken from Ref. [ll] and remeasured from 500 K down. As a result of thermal 
cycling/ annealing in this study (see text), the resistivity values are lowered by ~20 % from those corresponding values originally 
reported in Ref. Hll . 



Film 


2 flux 


p(300K) 


Pi 


Ei 


P2 


E 2 


P3 


E 3 


No. 


(SCCM) 


(fi cm) 


(fi cm) 


(meV) 


(fi cm) 


(meV) 


(fi cm) 


(meV) 


2 


0.02 


1.66 


0.079 


120 


1.89 


35.5 


2.63 


2.3 


3 


0.10 


2.72 


0.407 


113 


0.758 


39.2 


5.00 


3.1 


4 


0.15 


2.01 


0.189 


135 


1.30 


37.5 


5.88 


2.8 


5 


0.25 


5.89 


0.420 


135 


2.70 


26.3 


370 


5.7 


6 


0.50 


11.3 


0.588 


129 


4.79 


27.0 


333 


5.6 


7 


0.80 


26.0 


2.04 


121 


9.09 


37.7 


100 


4.3 



3 orders of magnitude lower than that in the single crys- 
tal. Such huge differences in resistivity provide a direct 
indication of the presence of far more numerous amounts 
of oxygen vacancies in polycrystalline films than in sin- 
gle crystal. (The shallow donor concentrations in the 
polycrystalline films will be estimated in the subsection 
HIB.3.) 

Our extracted E% values of «26-39 meV for the shal- 
low donors are in consistency with those previously re- 
ported for the ZnO materials! 16 ' 26 Our extracted E\ val- 
ues of «125±10 meV are in line with that (wllO meV) 
deduced by Tampo et ali^ in ZnO films grown by radical 
source molecular-beam epitaxy. However, these E\ val- 
ues are significantly lower than that («330 meV) found 
in the single crystal discussed in the above subsection. 
The reason why requires further investigation. One plau- 
sible explanation would be that, as compared with the 
single crystal, there are additional impurity levels which 
were unintentionally introduced during the growth of the 
polycrystalline films. On the other hand, since the Ei- 
conduction channel dominates the electrical-transport 
behavior up to ^400 K in our films, electrical measure- 
ments up to T sufficiently higher than 500 K would be 
highly desirable for an unambiguous determination of the 
values of the E\ activation energy and the associated 
charge conduction process.— 



VRH conduction processes previously observed below 
^90 K as mentioned in the Introduction, to be robust. 
This assertion is examined in the following. 
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0.25 0.30 0.35 0.40 0.45 

-1/4 -1/4 

T (K ) 

FIG. 4. (color online) Logarithm of resistivity as a function 
Q f j<-i/i £ or £ our ZnO polycrystalline films, as indicated. The 
symbols are the experimental data and the straight solid lines 
are least-squares fits to Eq. ©. 



2. Three-dimensional Mott variable-range-hopping 
conduction: 20-80 K 

As mentioned, our ZnO polycrystalline films had been 
measured previously in Ref. [Oj. In this work, we have 
extended the p(T) measurements on these films up to 500 
K. Notice that these measurements have led to thermal 
cycling/ annealing effect (in vacuum) on these films. As 
a result, their p(T) values are changed from the corre- 
sponding values originally reported in Ref. ITTI . In par- 
ticular, the p(T) values are largely reduced in those films 
grown under high oxygen atmospheres. For instance, the 
p(300 K) value decreased from 206 to 26 ft cm in the film 
No. 7. Nevertheless, we would expect the charge conduc- 
tion mechanisms, i.e., the Mott and the Efros-Shklovskii 



Figure |4] plots the variation of logp with 7 1-1 / 4 for four 
representative ZnO films, as indicated. The symbols are 
the experimental data and the straight solid lines are the 
least-squares fits to the Eq. ([3]) given below. Clearly, in 
every sample, there exists a linear regime from ~80 K 
down to ~25 K. This observation suggests that, in this 
T interval, the dominate electrical conduction process is 
due to the Mott VRH conduction mechanism in three 
dimensions i 2 ^ 

PM {T) = p Mo e^ 1/ \ (3) 

where pjv/ i s a temperature insensitive resistivity pa- 
rameter, and TM=18/[/csA r (i?F)C 3 ] i s a characteris- 
tic temperature. N(Ep) is the electronic density 
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TABLE III. Values of relevant parameters for six oxygen deficient ZnO polycrystalline films. pM Q and Tm are defined in Eq. ([3|. 
N(Ep) is the electronic density of states at the Fermi energy, ^?hop,Mott is the average hopping distance, £ is the Bohr radius of 
the shallow donors, and Whop, Mott is the average hopping energy. The values of Whop.Mott were calculated for a representative 
temperature of 40 K. 









N(E F ) 


-Rhop.Mott 




Film 




T M 


Whop, Mott 


No. 


(fi cm) 


(K) 


(J 1 m- 3 ) 




(meV) 



1.36 
0.61 
1.41 
1.15 
1.36 
4.18 



65 
3510 
854 
10400 
23600 
15200 



2.5xl0 48 
4.6xl0 46 
1.9 xlO 47 
1.6 xlO 48 
6.9xl0 45 
1.1 xlO 46 



2.89/T 1/4 
2.03/T 1/4 
3.79/T 1/4 
4.65/T 1/4 
4.16/T 1/4 



0.98 
2.64 
1.85 
3.35 
4.14 
3.68 



of states at the Fermi level, and £ is the local- 
ization length of the relevant electronic wave func- 
tion. The average hopping distance and the aver- 
age hopping energy of the conduction electrons are, 
respectively, given by J Rho P ,Mott=(3£/8)(T fl/ /T) 1 / 4 and 
Whop,Mott=(fcsT 1 /4)(rM/T) 1 / 4 .^ Our extracted values 
of Tm, and then calculated values of N(Ep), i?hop,Mott, 
and Whop,Mott are listed in Table III. In order to calcu- 
late these three quantities, we have used £«2 nm, the 
Bohr radius of the shallow donors in ZnO£ Our fitted 
-^-hop Mott values are significantly smaller than our film 
thickness (psI /an). Also, the criterion -Rhop,Mott/£ > 1 
for Eq. (J3j) to be valid is satisfied in the films Nos. 3 and 
5-7. On the other hand, this criterion is not satisfied for 
the films Nos. 2 and 4, because these two films are not 
resistive enough to lie sufficiently deep on the insulating 
side of the metal-insulator transition.— Furthermore, our 
extracted values of Wkop,Mott(40K)«l-4 meV are some- 
what smaller than the corresponding values reported in 
Ref. [ll|. This small decrease in Whop, Mott by ~35% from 
that in Ref. [ill can be readily understood as arising from 
the increased shallow donor concentrations in these films 
as a result of thermal cycling/annealing. Obviously, on 
average, a higher donor concentration will lead to a lower 
thermal activation energy involved in the VRH conduc- 
tion process. 



Our observation of the Mott VRH conduction in the 
T interval ^20-80 K in this work demonstrates that 
this carrier transport process in the ZnO material is 
generic and robust, while the values of the relevant pa- 
rameters may differ somewhat from run to run, because 
the amounts of donor concentration could have varied. 
At even lower measurement temperatures (<20 K), a 
crossover from the Mott VRH conduction process to the 
Efros-Shklovskii VRH conduction process 3 ^ should be ex- 
pected, as has previously been reported in Ref. [Tl|. Un- 
fortunately, the base temperature of our closed-cycle re- 
frigerator does not allow p(T) measurements down to suf- 
ficiently low T to illustrate this latter VRH conduction 
mechanism in the present work. 



3. Estimate of shallow donor concentration 

While the shallow donor concentration, ri£>, in the 
ZnO material is a very important quantity, it is dif- 
ficult to measure directly. However, we may estimate 
this quantity from the extracted p2 value which charac- 
terizes the ^-conduction channel in Eqs. (fT|) and {3J. 
At sufficiently high temperatures (in practice, above a 
few hundreds K), the shallow donors should have essen- 
tially all been ionized. These excited electrons would 
move in the conduction band and be in response to an 
externally applied electric field. The resulted electrical 
conductivity from this ^-conduction channel can ap- 
proximately be expressed by p^ 1 = {noe 2 T e ) / 'm* —nE>ep e , 
where T e is the electron elastic mean free time, m* is 
the effective electron mass, and p, e is the electron mobil- 
ity. Typically, the value of p e in ZnO at T above a few 
hundreds K can be extrapolated to a magnitude of ~10 
cm 2 /V s4 In the case of our as-grown single crystal, Table 
I lists a value p2 ~ 1x10 fl cm. Then, we may estimate 
'i_D=(p2e^ e ) _1 ^f 14 cm" 3 . In the case of our polycrys- 
talline films, Table II lists a value p2 ~ 1 Q cm. Then, 
we may estimate nD={p2Gp: e ) ^10 cm" 3 . These in- 
ferred shallow donor concentrations in the single crystal 
and in the polycrystalline films, respectively, are in good 
accord with the accepted corresponding values reported 
in the literature 4^ This degree of agreement provides 
a self-consistency check of our data analyses and the in- 
terpretation of the charge transport processes in ZnO in 
this section. 



IV. CONCLUSION AND SUMMARY 

We have investigated the electrical conduction proper- 
ties of ZnO single crystal and polycrystalline films in the 
wide temperature range 20-500 K. In the case of single 
crystal, we found that two types of thermal activation 
conduction processes dominate the carrier transport at 
T>110 K. As a result, a group of shallow donors and a 
group of intermediately deep donors, together with their 
individual ionization energies, have been inferred. In par- 
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ticular, we observed that the shallow donor concentra- 
tion is markedly affected by thermal cycling in vacuum, 
strongly suggesting the important role of oxygen vacan- 
cies as shallow donors in the ZnO material. In the case of 
oxygen deficient polycrystalline films, additional conduc- 
tion processes due to the nearest-neighbor-hopping con- 
duction and the Mott variable-range-hopping conduction 
mechanisms are observed. These two additional conduc- 
tion possesses originate from the existence of far more 
numerous amounts of shallow donors (oxygen vacancies) 
in polycrystal films, which led to the formation of a par- 
tially filled impurity band. Taken together with our pre- 
vious work (Ref. Illl ). this present study provides a fairly 
complete picture of the overall electrical conduction pro- 
cesses in the ZnO materials. 

A summary of the overall electrical conduction pro- 
cesses in the ZnO materials: We would like to give a 
brief summary of the rich electrical conduction processes 
in the intensively studied, scientifically and technologi- 
cally alluring ZnO materials. As the temperature pro- 
gressively increases from liquid-hilum temperatures up 
to 500 K, one expects to see, one process by one process, 
(1) the Efros-Shklovskii variable-range hopping conduc- 
tion, (2) the Mott variable-range hopping conduction, (3) 
the nearest-neighbor-hopping conduction, (4) the ther- 
mal activation conduction from the shallow donor levels, 



and finally (5) the thermal activation conduction from 
the intermediately deep donors levels. The conduction 
processes (l)-(3) take place in the shallow donor impu- 
rity band, while the conduction processes (4) and (5) take 
place in the conduction band. These conduction mech- 
anisms are generic and robust in unintentionally doped 
n-typc ZnO materials which lie on the insulating side of 
the metal-insulator transition^ However, the tempera- 
ture interval within which each conduction mechanism is 
found and the extracted value of the associated parame- 
ter may vary more or less from experiment to experiment. 
Such variations can readily arise from the fact that dif- 
ferent experiments usually involve differing amounts of 
donors and differing donor level distributions in the sam- 
ples. 
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